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(54) Method and apparatus for spectral analysis 

(57) A method of analyzing an optical image of a 
scene to determine the spectral intensity of each pixel of 
the scene, which includes collecting incident light from 
the scene; (b) passing the light through an interferom- 
eter which outputs modulated light corresponding to a 
predetermined set of linear combinations of the spectral 
intensity of the light emitted from each pixel; focusing 
the light outputted from the interferometer on a detector 
array; and processing the output of the detector array to 



determine the spectral intensity of each pixel thereof. If 
the interferometer is of the moving type scanning in one 
dimension is required where the detector array is one 
dimensional, and no scanning when the detector array 
is two-dimensional. If the interferometer is of the non- 
moving type scanning is required in one dimension 
when the detector array is two-dimensional, and in two 
dimensions when the detector array is one-dimensional. 
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Description 

FIFI n AND BACKGROUND TMF INVENTION 

™. p,e Se n, ^ — »a - = -^Sr*S^C 

lyzing an optical .mage of a scene to det f™* ^J D ™ 'J^nert wavelengths, and detect the spectrum. An 

is covered by a one-dimensional scanner and the remaining dimension of the j ar ny wnen a o 

fie.d of view is scanned mechanically in two - S ^ enures S 

time only to sample the different wavelengths of a angle port In ^ J J it would be 

each detector sees only the contribution of a single pixel at a single wavelength at any i.me. o 
impossible to separate the spectra of each pixel. rii.advarrtaae that most of the pixels of one 

However, conventional slit-type imaging ^f° m *^ 
frame are not measured at any given tone, even though the tore ophes o ' * e * e ™ , ; e J avai|able information since, 
from a., of them simultaneously. Thus, the «™« rt r^ Jven "me and reaching a par- 

except for one wavelength^ most oMh. ^ZTJm^I is required to obtain the nec- 

S~ a tfon S-iS^ - £ sign y a,-to 9 -noise ratio (sensitivity) is substantial, reduced for 

3 91 An o^cT^rpreTent invention is to provide a nove. method and apparatus for spectra, analysis of images which 

have advantages in the above respects. ^ athoH anc i annaratus for spectral analysis of images 

the required frame time and/or to substant.ally increase the signal-to-no.se ratio, as comparea 
type of imaging spectrometer. 

SI IMMARY OF THF INVENTION 

interferometer, as well as of the non-mov.ng type wherein i the OPD. is var.eo I w.tn ne g ^ 
radiat.on. Thus, in the moving type interferometer, at each instant each J^^J^™^ time . When the 
its signal is a linear combination of the spectral content of the ra ^ ^^relevantLear combinations 

been scanned at all relevant linear combinations of the spectral ^"J- h Fabry-Perot and Michel- 

For purposes of illustration, the invention is described below as Michelson and 

son interferometers as examples of the moving-type interferometers, and as .mplemented by use 

Sagnac interferometers as with the above method. 

The inventor, also provides s ^ t "* * r j*^ differ from the conventional slit-type imaging 

The methods and apparatus in accordance with the ^ n jtjn tne ^ 

« spectrometer by utilizing an interferometer as de ^' b ^ b ° ve - '"^^ Such methods 

lected energy with an aperture or slit, thereby ^ stant f *' n <^ analyzed, thereby 
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trial and Planetary Remote Sensing, SPIE Proceedings, Vol. 750. p. 140 (1987). 

Let n be the number of detectors in the linear array, m x m the number of pixels in a frame and T the frame time. 
The total time spent on each pixel in one frame summed over all the detectors of the array is: 



nT/m 2 



By using the same size array and the same frame rate in a method according to the present invention, the total time 
spent summed over all the detectors on a particular pixel is the same. nT/m 2 . However, whereas in the conventional 
grating or prism method the energy seen by every detector at any time is of the order of 1/n of the total because the 
wavelength resolution is 1/n of the range, in a method according to the present invention the energy is of the order of 
unity because the modulating function is a sinusoidal (Michelson) or similar periodic function (low finesse Airy function 
with Fabry-Perot) whose average over many periods is 50%. Based on the standard treatment of the Jacquinot advan- 
tage (or multiplex advantage) described in interferometry textbooks, it is possible to show that devices according to the 
present invention have measurement signal-to-noise ratios which are improved by a factor of n in the infrared range 
(background limited performance) and by the square root of the ratio of the signal at a particular wavelength to the aver- 
age signal in the spectral range, at wavelengths of a narrow peak in the visible range (signal photon noise limited per- 
formance). For a mathematical treatment and definition of the Fabry-Perot interferometer and a definition of the Airy 
Function, see Max Born and Emil Wolf. Principles of Optics, Pergamon Press. 1980. p. 329. 

In all the embodiments of the invention described below, all the required optical phase differences are scanned 
simultaneously with the spatial scanning of the field of view in order to obtain all the information required to reconstruct 
the spectrum, so that the spectral information is collected simultaneously with the imaging information. 

The invention can be used with many different optical configurations, such as a telescope for remote sensing, a 
microscope for laboratory analysis, fiber optics for industrial monitoring, and others. In addition, any wavelength range 
can be selected with appropriate filters and optics. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The invention is herein described, by way of example only, with reference to the accompanying drawings, wherein: 

FIG. 1 illustrates a conventional (prior art) slit-type imaging spectrometer; 

FIG. 2 is a block diagram illustrating the main components of an imaging spectrometer constructed in accordance 
with the present invention; 

FIG. 3 is a diagram illustrating an imaging spectrometer constructed in accordance with the present invention 
based on the use of a moving-type interferometer, namely, a Fabry-Perot interferometer with scanned thickness; 

FIG. 4 illustrates a Michelson interferometer with retroreflectors used as a non-moving type interferometer in the 
imaging spectrometer of Fig. 2; 

FIG. 5 is a diagram illustrating another moving-type interferometer, namely, a Michelson type interferometer, which 
may be used as the interferometer in the imaging spectrometer according to the present invention; 

FIG. 6 illustrates another non-moving type interferometer, namely, a modified Sagnac interferometer, as used in an 
imaging spectrometer in accordance with the invention; 

FIG. 6A illustrates yet another non-moving type interferometer, namely a modified Sagnac interferometer, as used 
in an imaging spectrometer in accordance with the invention; 

F!G. 7 schematically demonstrates that the rotation of a ray at the entrance to an interferometer such as that of Fig- 
ures 6 and 6A also rotates the ray at the output by the same angle and changes the OPD between the correspond- 
ing split beams; 

FIG. 8 schematically demonstrates that the same effect as demonstrated in Figure 7 can be achieved through the 
rotation of the interferometer; 

FIG. 9 illustrates the invention as embodied in a modified Michelson imaging spectrometer, as an example of 
another moving type interferometer, with focusing on the mirrors through microscope objectives. 
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FIG. 10 shows a modified Sagnac interferometer as in Figures 6 and 6A but where the scanning is effected by a 
movable element of the interferometer; 

FIG. 1 1 shows an interferometer as in Figures 6 and 6A, but with an additional optical plate made of a light trans- 
mitting material at 90° to the beamsplitter; 

FIG. 12 shows an interferometer similar to that of Figure 1 1 but with a different orientation of the various compo- 
nents, the additional optical plate 138 is also at 90° to the beamsplitter; 

FIG. 13 shows a Michelson-type interferometer similar to that shown in Figure 5 but with a two-dimensional detec- 
tor array; 

FIG. U shows a Sagnac interferometer as in Figure 1 but with a scanning mirror rather than a beamsplitter; 
FIG. 15 is similar to Figure 14 except that the scanning mirror is single-sided. 
DESCRIPTION OF THE PREFERRED E MBODIMENTS 

For purposes of better understanding the present invention, as illustrated in Figures 2-1 5 of the drawings reference 
is first made to the construction and operation of a conventional (i.e.. prior art) slit-type imaging spectrometer utilizing 
a two-dimensional array of detectors as illustrated in Figure 1 . cuetom 

Thus, the prior art slit-type imaging spectrometer as illustrated in Figure 1 comprises a collector i opt cal systenv 
such as a telescope as indicated at 2. for collecting the incident light from the scene, schematjcaHy mdica ted I at 4 and 
focusing the substantially parallel light of the scene onto a first focal plane occupied by a slit 6 to def .ne J^*^ 
The light exiting from slit 6 is collimated in a collimator lens 8 and is passed through a transm.ss.or, or reaction grating 
1 0 to separate the various wavelengths. The output from grating 1 0 is focused by a focusing lens 12 onto a two-d.men- 
sional detector array 14 in a second focal plane. The output of detector array 14 is fed to a s.gnal processor i* 

in the two-dimensional array of detectors 14 illustrated in the prior art imaging spectrometer of Figure 1 , the move- 
ment of the system (e.g., when embodied in an aircraft) effects the scanning along one dimension. The scanning along 
the second dimension is effected by the slit 6 which is oriented perpendicularly to the direction of movement of the sys- 
tem. The slit 6 thus assures that each detector within the array 14 sees only the contribut.on of one pixel at a s.ngle 
wavelength at any time; this is necessary to separate the spectra of each pixel. 

As ment.oned earlier, the disadvantage of the prior art method illustrated in Figure 1 is that most of the pixels of one 
frame are not measured at any given time even though the telescope (2 or other collecting optos) actua y collecte 
energy from all of them simultaneously. As a result, the required frame time is significantly ■""eased, and/or ^ signal- 
to-noise ratio (sensitivity) is substantially decreased with respect to a system, if it existed. wh.ch d.d not have the need 

,0f ^rel is a block diagram illustrating the main components of an imaging spectrometer constructed in accordance 

with the present invention. , n . , nna 

Thus, the imaging spectrometer of Figure 2 includes: a collection optical system, generally designated 20, a one 
dimensional or two-dimensional scanner, as indicated by block 22; an optical path difference (OPO) generato. -or inter- 
ferometer, as indicated by block 24; a one-dimensional or two-dimensional detector array, as indicated by block 26, and 
a signal processor and display, as indicated by block 28. • . ^„„ r , K 

Throughout it is intended that the radiation, such as light, to be analyzed can come from a w.de variety of sou ces. 
For example, the source may emit radiation spontaneously or reflect or transmit radiation from a lamp or other i urn,- 
nated object. In addition, with proper illumination, such as UV of laser, and with proper means of preventing the illumi- 
nating wavelengths from reaching the spectral imager, fluorescence or Raman spectral imaging measurements can be 
performed, in order to obtain different information about the about the object or objects in question in i each i case. 

A critical element in the novel system is the optical path difference generator or interferometer 24 .. wh.ch outputs 
modulated light corresponding to a predetermined set of linear combinations of the spectral intensity of the ight emitted 
from each pixel of the scene to be analyzed. The output of the interferometer is focused onto the detec to r ar ay 26_ 
Thus all the required optical phase differences are scanned simultaneously with the spatial scanning of the field of 
view 'in order to obtain all the information required to reconstruct the spectrum. The spectrum of all the p.xels in the 
scene is thus collected simultaneously with the imaging information, thereby permitting analysis of the image in a reai- 

time A m methad and apparatus according to the present invention may be practiced in a larg. ^ ot 
Specifically, the interferometer used may be of either the moving or the non-movmg type and the detector array , ma* 
•ndependently of the type of interferometer, be one- or two-dimensional. When the interferometer is of the rnov.ng type 
and the detector array is two-dimensional, no scanning is required, except for movement of the .nterferometer which .s 
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an OPD scan When the interferometer is of the moving type and the detector array is one-dimensional, spatial scan- 
ning in one dimension is required. When the interferometer is of the non-moving type and the detector array .s two- 
dimensional. OPD scanning in one dimension is required. When the interferometer is of the non-mov.ng type and the 
deTertor array is one-dimensional, scanning in two dimensions is required, with one dimens.on relating to a spatial scan 
; while the other relates to an OPD scan. 

Figure 3 illustrates one form of imaging spectrometer constructed in accordance with the, ™« 
spectrometer is based on the use of a moving type interferometer in which the OPD ,s vaned to modulate the hght. 
namely a Fabry-Perot interferometer 33 with scanned thickness. 

ThuJ the imaging spectrometer illustrated in Figure 3 comprises a source or scene, generally des.gnated 3CUo be 
,o analyzito dSineme spectra, intensity of each pixel of the scene. The scene 30 may be a source o incoherent 
nofmono^hromatic radiation. K may be at a distance for remote-sensing applications, in which case the collection opti- 
ca. system schematically indicat ed at 31. would be a telescope; alternatively. *• scene 30 may be close for micro- 
scoDic ana ysis. in which case the collection optical system 31 would be a microscope. 

Optical system 31 provides an output to a one-dimensional mechanical scanner, e.g.. a mirror scanner as md-cated 
,5 at 32 to scan the scene. The output from scanner 32 is fed to a Fabry-Perot interferometer 33 havmg an etalon made 
of too plane, parallel reflectors spaced at a distance V from each other. In this example, the spaang distance d ,s 
variable by using a mechanical scanner, in this case a piezoelectric scanner 34. nn on 

The output from the Fabry-Perot interferometer 33 is fed through a refocusmg optical system 35 onto a one^imen- 
sional array of detectors 36 whose outputs are fed to a signal processor 37. If a two-dimens.onal detector array .s used. 

20 OptiS fystemlreran afoca. telescope or microscope, produces a substantially paral.el beam (i.e., exactly 
oarallel or having a very large F/No) at its output, because in this way every detector within the array 36 corresponds to 
a sTngt optica, phase difference through the etalon 34 of the Fabry-Perot interferometer 33. The optica, system ^31 can 
be eTerTefraive or ref.ective. The etalon 33 of the interferometer is at 90" to the optica, axis of the system, .t w... be 

oc noted that no use is made of a field-of-view limiting aperture or slit 

Ihefmaging spectrometer i..ustrated in Figure 3 operates as follows: A detector element i of the array 36 rece.ves 
radiation from a beam which traverses the etalon 33 at a fixed angle ( W between its line of sight and the norma, to the 
etalon. and therefore it sees this radiation always through an optical phase difference b„ given by: 

6 i= 2 rt 2d(n 2 -sin 2 *,) 0 - 5 /*. (1) 

where >. is the wavelength of the radiation considered, and n is the index of refraction of the air between the plates^ 

The total radiation teaching detector i k at any given time from a specific pixel k of the radiation ,nput 30 s the con- 
volution <rf that pixel's spectrum with the Airy function, as follows (Max Born and Emi. Wolf. Prmcp.es of Optics, Perga- 
:~5 mon Press, 1980. page 327): 

l ik = J {!„<*) /(1 +Fsin 2 [8 i (X)/2])}dX (2) 
where the integration is between X, and X 2 . and 

X, , >. 2 = lirni ts of the spectral range 

\ k n ) = spectral intensity of the source at pixel k 

F = coefficient of finesse, is related to the reflectivity of the etalon R. by: 

F = 4R / (1 - Pi) 2 

i k = that particular detector which images pixel k through an optical phase difference 6,{X) 

Following is one of many ways of scanning the field of view and the thickness of the etalon 33. 
Suppose the array 36 is composed of a linear set of N detectors, whose signals can be mon-tored *muKaneous y 
and fndependently. Suppose M performs a raster type scan of m lines (larger than N) and that the p ane of the paper n 
Figure 3^s the vertical direction. Every time M has scanned one horizontal line, the thickness d I of the «-£n33 J 
incremented by the piezoelectric scanner 37 one step in synchronization with the vertical scanner 32 starting from d-0 
unt *N lines are scanned, and N thickness steps are made. At this moment the etalon th,ckness ^£ 
inal value, and the thickness scanned again. The process is repeated until the scanner has scanned one complete 

''^Except for a marginal region of N pixels at the top and at the bottom of the field of view, all the pixels of the f,e.d of 
view are measured with N optical phase differences by different detectors. All the ^tector signals a 
recorded at a high rate, such that all the needed information is collected and fed to the s.gnal processor 37 to recon 



45 



5 



10 



15 



struct both the image and the spectra of all the pixels. ,. . rtre , . 

As mentioned above, another possible conf iguration is one which includes a two-d.mens.onal array of detectors. ,n 
thi«i case the same idea applies, but one spatial scan is saved. 

For eTarile tf £ alay is an N x m matrix in a "push broom" approach similar to that descnbed Ion Page 42 of 
John £ We.Sn, imaging Spectrometers for Terrestria. and Planetary Remote Sensing. SP.E Proceed.ngs. Vol. 750. 
o 1 An MQft-n the N lines and m columns correspond to the same matrix in object space. 

" his SnfTauration the Fa^ry Perot thickness or optica, phase difference is kept fixed for the t.me o^egrat-on of 

One divides the spectral range of interest X, to X 2 . into N intervals. If one approbates the .ntegral of Equafon (2) 
as a sum over N wavelength intervals. Equation (2) is the product of an N x N dimens,ona. matnx g,ven by. 

A.,-1/11 +Fsin(26J2)] (3) 



and an N dimensional vector l^), where j scans the N wavelength intervals, and * * ■ {** ^ J^jjE 

If one inverts the matrix A. and multiply it by the vector l jk (i k =1 .... to n). one obtams the vector wh,ch is the spec 



20 If one inverts the matrix A 

"^wlth're'wect to spectral resolution, consider two quasi monochromatic sources placed at the same pixel. 

narrol peaks of the Airy function. This is a desirable station, since in any case, because of the ,mag,ng, a h.gh reso 

of the Born and Wolf pubiication s,< ho.ds. because if one sets, 
30 sin(c/4) = ( c /4) - 1/(F) 0 5 = 1/0.3 - 0.3 ( 4 > 

one sees that one makes an error of about 10% in approximating sin e/4 with eM. 
So one has in this case, 

55 

^=7i/2 (F) 0 J = 5 (5) 

4C and the resolving power is, 

X/AX ~ 2 ^nd/ * ( 6 ) 



In order to 



get an idea of the order of magnitude, and check consistency, let it be assumed that, 



N = 50 detectors in the array (7 > 
=3u/50 = 0.06m (8) 



and. 

k = 2 to 5u (9) 

Therefore, d is the range, 
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d=X 2 /AX(2ir n )- (10) 



Taking n=l for air, Equation (10) gives, 

d = 10/(0.06 x 10) = 16^ 

The range of d is such that it scans the same range of optical phase differences as the wavelength range; therefore. 

d 2 a 1 -d 1 /x 1 -dAA., -d/x 2 < 12) 

or 

da-d, -d[1 -X^X,,]* 16x0.3 = 5^ ( 13 > 

Therefore the thickness steps are of the order of. 

S^i/50 = 0.1 p. (14) 



35 



In the 8 to 1 4^ spectral range one has. 

d = 11 2 / (0.12x10) = 100m (16) 



25 AX = 6/50 =0.12}i (15) 



and, 

d 2 -6, = 100[1 -8/14] = 43p < 17) 
Therefore, the steps of thickness are of the order of 1 ^ In the visible range 0.4 to 0.8* AX = 0.008 M , and therefore. 

d = 0.6 2 / (0.008 x 10) = 4.5u ( 18) 



From (18) 



and therefore the steps in d are 



d 2 -d, =4.5[1 -1/2] = 2.25ji 



(19) 



45 



2/50 = 0.04 M (20) 



55 



in summary the significant features in the system illustrated in Figure 3 include: i) the special matching of the inter- 
JSii and f inesse, with the detector array size and number ^.^^^J^^;. 
and ii) the synchronization between the thickness scanning and the spatial scannmg. to obtain the spectral and the spa 
tial information simultaneously in the time one frame is built. , ra rmn* the 

Many other scanning configurations may exist, which do not require synchronization. ^^JT^^J^ 

SC ° P ,^I b 6 e sTwh^S the movng type interferometer as illustrated , Figure 

eter at a small angle (*-0) to the optical axis undergoes an optical path difference which va nes as * or J'O^POwe 
OH- All the spectra, information in all the pixe.s may be collected by scanning the ^^^^JJ^Si 
zation with the scene scan, at the end of which, every pixe. has been measured thmughjJM 
by different detectors. By careful bookkeeping and by applying the appropriate matrix mvers.on (such as Four.er 
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formation), the spectrum of every pixel may be calculated. The bookkeeping is ^^^^^Z 
gather the information of different OPD s of one pixel at different times Thus in the time of f O msec the usual frame 
me of a standard video), a spectrum may be measured for every pixel of a standard v,deo frame. Th.s .s of the order 
of 100 resolution points per pixel, with a typical matrix of 500 x 500 pixels per frame. 

in th» instruction of Fiaure 4 instead of using a moving-type interferometer wherein the OPD (optical path differ- 
ence? varied S ^Sement inlhelnterferometer (namely, changing the spacing d between tt^J^ 
the Flbry^eSt^ometer illustrated in Figure 3). the invention may also be implemented by using ^moving 
tne ra ° r y re ™ 1 m opD j s wjth the angle 0 f incidence of the incoming radiation. Figure 4 illus- 

S IT 2. «» . « ^ a ». «■ 

"'"Thus. M intele-cme.* «uMM in Figure 4 indudss a beamsplitter '"-7^* 

^ ^S'i^he Micheison type interferometer illustrated in Figure 4. the beam which is paraile. to the optical axis of 

„rJ in Fiaure 3 each element of the array will receive light which underwent different OPD s between the two arms. 

to the beamsplitter 40. and compensator 42 perpendicular to it. the optica, path 
difference between the two arms can be shown to satisfy for small f s. the following relationsh.p: 

OPD(<t>) = 2d sin <t> 0 [1- cos<j> 0 /(n 2 -sin 2 <t> 0 ) 1 <t> {21) 

with d, n - angle between the beamsplitter 40 and the optical axis. For simplicity can be taken as 45-. but this is not 
Tssential nTs thetcS of refraction of the two compensators. 41 and 43. d is the compensator thickness, and ♦ is the 

an9 for e^ "Cir^nge in OPD between two adjacent detectors is -1* and the change in . directions 
be tween ^ l heTe.ds of Ve^of two neighboring detectors is - 1 mrad, d turns out to be of the order of 1 mm. In fact. 

d = 1n / {2(2) 05 /2 [1 - ((2)°' 5 /2)/(4-l/2) 05 ) 10" 3 ] = 1 mm 

Note that if compensator 41 is removed and beamsplitter 40 is half-reflecting on the n ^'e surface on the side of 
the incoming beam, then one obtains another viable configuration with a compensa ed optica I ax * 

With respect to the spectral resolution of the Figure 4 configuration, assume that the spectral range of the rad.at.on 
reaching the detector is limited to. 
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hv . ,a\na a suitable filter or because of the transmission properties of the optics. non 
^ "BeLuse S the Nyquist samp.ing theorem, in order to avoid aliasing, the i"*"^^ 1 ^^^^ 
steps which are no larger than X,I2. For the purpose of illustration, one can assume that the *J^**™%^£ 
Spd subtended by two adjacent detectors. Since one period is an OPD change of one wavelength, the max.mum d.f 
ference in OPD seen by two adjacent detector elements must be, 

Assume now that the incident radiation consists of two lines at wave numbers v and v + Av with an intensity of l 0 - 

The resultant intensity for an OPD of x is: 

l(x)=0.5l 0 (1 +2cos2;rvx)+0.5l 0 (1 +2cos2it(v+ Av))x) (23) 
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If one subtracts the constant term, one obtains, 

l(x) = l(x)-l 0 < 24 > 
= 1 0 [cos(2jtvx) + cosit(v+Av)x] 
= 21 o [2:t(v+Av) x] cos Tt AVX 

Similar to the Rayleigh criterion, one defines the two line "resolved" if the ratio, 

i(x)/l(0) 

is lower than a predetermined amount, say 0.9. This gives a condition for the maximum value of x, x max . needed to get 
the defined resolution. 

The value of x max which satisfies the above condition is: 

cosrr Avx max *0.9 ( 25 ) 



or, 

20 x maxS 0.143Mv (26) 

Now let N d be the total number of the detectors in the array. Then: 

x m ax^N d X,/2 (27) 
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and, 

N d X,/2 =0.143/Av (28) 



or. 



Av = 0.286/(X 2 N d ) (29) 

For vg = 



Av/v, = Av k 2 = (0.286/N d (>. 1 /X. 2 ) 



As a numerical example: 

if = 5um. k 2 = 2um and N d = 100 



Av/v 2 = 0.7% 



In summary, the uniqueness of the system illustrated in Figure 4 is represented by: i) the OPD ,s a linear function 
of the angle of incidence of the incoming radiation on the interferometer, so that different detectors of the array see it at 
different OPD s This fact, combined with the spatial scanning and proper bookkeeping, allows the .nterferograrn or Fou- 
rier transform of the spectrum of each pixel to be measured simultaneously with the image information. Again. *e scan- 
ning in the case of a non-moving interferometer can be two-dimensional or one-dimensional, depend.ng on whether the 
detector array is one-dimensional or two-dimensional, respectively. -■■„»„ 

Figure 5 illustrates an imaging spectrometer including a Michelson interferometer but of the mov.ng type, ami ar to 
that of Figure 3 namely, wherein the OPD varies with moving an element of the interferometer. In the spectrometer of 
Figure 5 the light from source 50 is collected by the optical collection system 51 and is collimated onto a scanning mir- 
ror 52 before being passed through the beamsplitter 53 which splits the beam into the two arms. One arm includes a 
compensator 54 and a mirror 55 and the second arm includes merely a mirror 56. The light from the two mirrors ,s 
passed via beamsplitter 53 and a focusing lens 57 onto an array of detectors 58. „„„ ir ,n « our 

In the spectral imager illustrated in Figure 5. if the detector array is linear then the d.rect.on of the scann.nc , is per- 
pendicular to the direction of the linear array and the scanning mirror 52 scans the scene in one d.mens.orv Sea inner K 
is not needed if the detector array is two-dimensional. Scanner 59a controls the distance between mirror 56 ^ beam- 
splitter 53 By scanning mirror 56 in the way of the traditional Michelson interferometer, while mirror 55 is st at jonary: ne 
OPD of the two arms is varied simultaneously for all pixels in the scene. The compensator 54 ensures that the central 
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rSS=3=35=3E3S£===S5 

(0 - 60° in Figure 6 and 45" in Figure 6A) goes through the .nterferometer w.th an OPD - 0. whereas ray 
incident at a general angle P - 6 undergoes an OPD given by the following: 

OPD( P .e.t.n) = t[(n 2 -sin 2 ( P+ e)) 05 - (n 2 -Sin 2 ( P -e)) 05 + 2sin P sin9] (30) 

where, 

; p = is the angle of incidence of the ray on the beamsplitter ((60- in Figure 6 and 45° in Figure 6A) 

angular distance of a ray from the optica, axis or Sagnac rotation angle with respect to the centra, posit.on 
t = the thickness of the beamsplitter 
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the index of refraction of the beamsplitter 
Equation (30) may be approximated by the linear equation: 

OPD(p.e.t.n) = 2tsinp [1-(n 2 -sin2p)' 05 cosp) 9 (31) 
What follows is valid also for the conf iguration of Figure 4 and Equation (21 ). 

It follows from Equation (30) that by scanning both positive and negative angles with respect to the central pos.t.on. 
one can get a double-sided interferogram for every pixel, which helps eliminate phase errors giving more accurate 
°esuS?n me Fourier Transform calculation. The scanning amplitude determines the maximum OPC .reached, wh.ch is 

efatS to he spiral resolution of the measurement. The size of the angular steps determines the OPD step which*. 
In turn Sctated by the shortest wavelength to which the system is sensitive. In fact, according to the Nyqu,st theorem, 
this OPD step must be smaller than half the shortest wavelength to which the system .s sensitive. , . _ . 

Another parameter which should be taken into account is the f into size of a detector element in l^^^SE 
the focusing optics, the element subtends a finite OPD in the Sagnac interferometer which has the effect of convolving 

Ee Xerogram with a rectangular function. This brings about, as a consequence, a reduchon of system sens,tv ty at 
sho t wavelengths, which drops to zero for wavelengths equal to or below the OPD subtended by the e enient. For thjs 
reason one must ensure that the Modulation Transfer Function (MTF) condition is satisfied. ,e that the OPD sub- 
tend by a detector element in the interferometer must be smaller than the shortest wavelength at wh.ch the instru- 

men HQure n 9 t l e strates an addrtional imaging spectrometer including a Michelson interferometer of the moving type. 
simi.aMo tha Lstrated in Figure 5, but constructed so that the light is focused on the reflectors through microscope 
Sves Thus the spectrometer illustrated in Figure 9 includes a light source 71 and a m,croscope opt-ca system 72 
SST. condenseHens 72a, for focusing on the object 72b. Optica, system 72 further includes » 
for directing the light through a scanner 73. through a beamsplitter 74 to the two reflectors, 75 and 76. and then through 
riSS?^ " o^o a Ltector array 78. The detector array 78 in configuration of Figure 9, as well as in the previ- 
ously-described configurations, may be a CCD (charge coupled device). „ incnntnooriPd 
The scanner 73 scans the field of view in one dimension if the detector array is one-d.mens.onal. It is not needed 
at all if the detector array is two-dimensional. Reflector 75 and the microscope system move together along the optical 
axis to pTovfde scanning of the OPD in synchronization with the scanner 73 (if present), as descr.bed above w.th 

respect to Figure 5. 

As described above, the OPD for every pixel can be scanned m two ways: 

(1) A fixed type interferometer is used, in which the interferometer is kept in a fixed position and the scene , is 
scanned with an optical scanner. Since each detector in the detector array sits in a pos,t,on on the focal P ane from 
Sen it sees the radiation through a fixed OPD. which OPD varies from de tector ^^^T^^ 
ensures that each pixel is seen by all the detectors, and therefore through all the OPD s at the end of the scan. At 
the end of the scan, the Fourier Transformation can be applied to all the pixels; and 

(2) A moving type interferometer is used, wherein the scene is kept fixed, and an element .n the interferometer o 
hi who°e interferometer itself is scanned in such a way that now each detector always sees the , same pixel, bu 
hfough different OPDs at different times. At the end of the scan, all the pixels w,l. have been 

the OPD s. The Fourier Transformation can therefore be applied to the signal function measured by every detector. 

Where the detector is a linear array, a one-dimensional optical scan of the scene must be added to the interferom- 
eter ^an. if he entire scene is to be analyzed. In the case of a line scanner (where the scene » ^ 
area) no additional optical scan is needed, except for the interferometer, but there are many work.ng conf.gurat.ons 
depending on the dimensionality of the detector array and the type of interferometer used. 

The decsion as to which method is preferable should be based in each case or , the ^^^JS^ a 
metric limitations For example, in method (1). the spectral resolution is lim.ted by the number of detectors and bythe 
Sta OPrranTe allowed by'he geometry of the setup. In method (2). depending on the ^^^^^ 
is possible to reach very large OPD s. Since the spectral resolution atta.nable .s h.gher when a h.gher max.mum OPD 
ic scanned it is clear that the second method is more suitable when high resolution is required. 

^he^cSiment of Figures 10-15 are based on the second method, namely, the one wherein the scene .s fixed 
and the OPD scanning is effected by a movable element of the interferometer or the .nterferometer itself, and on a th rd 
met od wh^hTs bas'ed on a fixedinterferometer with an additiona, compensating ^ 

sided), in order to keep the scene stationary, and scan the interference fringes across the scene jn the s J°^ a ^ third 
methods, the light is passed through a moving type interferometer, in which the OPD JJJ 
moving an element o* the interferometer or a compensating m.rror. or the ^^^^^^ 
each detector sees a linear combination (Fourier transform) of the spectral content of the light em.tted from each pixel. 
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and it looks at a different pixel of the scene, and that when the scanner completes one scan, the complete scene will 
have been scanned at all relevant linear combinations of the spectral content. 

In all these methods the scene must still be scanned in one dimension if the detector array is linear. 
In the drawings Figures 10-15 illustrate further implementations of the second method and two .mplementat.ons of 
the third method, i.e. wherein the scene is kept fixed and the interference fringes scan the scene. 

Figure 10 illustrates an imaging spectrometer of a type similar to the modified Sagnac interferometer .llustrated m 
Figure 6 except that it is based on method (2) described above, namely, effecting the scanning by a movable element 
of the interferometer, rather than method (1) above, wherein scanning is effected by scann.ng the .mage. 

Thus the imaging spectrometer illustrated in Figure 10 is similar to that illustrated in F.gure 6. in that rt .ncludes an 
ootical image source 110 in which all the spectral information in all its pixels, after being collimated by the opt.cal col- 
Son system 1 1 1 . is scanned by a mechanica. scanner 1 1 2. This scanner is needed only if the detector array ,s hnear 
Ois two-dimensional, it is not needed. The light is then passed through a beamsplitter 1 1 3 to a f.rst reflector 114. and 
then to a second reflector 115, which reflects the light back through the beamsplitter 113 arxl then through a focusing 
lens 116 to a detector array 117. In the apparatus illustrated in Figure 10. however, the beamsplrtter 113 is rotated abou 
an axis perpendicular to the optical path (and perpendicular to the paper) to effect the scanning. Thus, th.s rotat.on o 
the beamsplitter changes the OPD through which a collimated beam from a pixel reaches a specrf.c detector: without 
changing Vne detector on which it is focused. It is thus clear that a scan of the beamsplitter angle with respect to the 
oDtical axis causes every pixel to be scanned through an OPO range. 

Thus in method (2) above (scanning being effected by the interferometer), when one looks at the focal plane where 
the detector is placed, one can see that the image of the scene is superimposed on a number of moving '^rfererice 
frinqes or alternates between dark and bright scene, depending on the type of interferometer. In method (1) above ( he 
scene scanning method), one can see that, as the optical scanner moves, the image of the scene scans through the 
focal plane, while the fringes remain stationary on this plane. Therefore, each point on the scene goes through positions 
on the focal plane, which correspond to different OPD's. At the end of one scan, all the needed information « collected 
lor each pixel. The results of the two methods are therefore similar, but the arrangements based on the second method, 
provide a number of advantages as described above. 

Figure 11 illustrates an imaging spectrometer similar to that of Figure 10 above or to F.gure 6. except that it, ncludes 
an optical plate 128 of a light-transmitting material, which increases the maximum OPD and thereby '"creases he 
spectral resolution of the apparatus. In this case, the optical plate 128 is located between the beams,: J ltter 123 a " d ^ 
second reflector 125, which reflects the light from the first reflector 124 back through the beamsphtter 123 and then 
through a focusing lens 126 to the detector array 127. AMararti „ riori t a 
Figure 12 illustrates an imaging spectrometer very similar to that of F.gure 2. except for a sl.ghtly different or.enta- 
tion of the first reflector 1 34, the second reflector 1 35, the optical plate 1 38. the beamsplitter 1 33, and the detector array 
137 including its focusing lens 136. ■ „ = 
Figure 13 illustrates an imaging spectrometer including an interferometer of the Michelson type, similar to F.gure 5. 
in this case, however, a two-dimensional detector array 147 is used, and therefore, a one-dimensional scanning .s suf- 
ficient rather than a two-dimensional scanning as in the apparatus illustrated in Figure 5. Thus, in Figure 13, the scan- 
ning is effected by mirror 146 controlled by scanner 149a, to change the distance between the mirror and the 
beamsplitter 143. and thereby to vary the OPD. Thus, the system illustrated in Figure 13 obviates the need for the syn- 
chronization circuit 59 and scanning mirror 52 included in Figure 5. 

Fiqures 14 and 15 illustrate two further embodiments of the invention which provide advantages over the embodi- 
ment illustrated in Figure 10. Thus, in the apparatus illustrated in Figure 10, based on the Sagnac interferometer with 
rotating beamsplitter, the distance between fringes varies as a function of beamsplitter angle. Th.s makes the transfor- 
mation algorithm inconvenient and difficult to use. In the embodiments illustrated in Figures 14 and 15. however, the 
fringes scan the image without changing the distance between them. In this way. a constant fringe d.stance ,s attained, 
with the possibility of reaching higher spectral resolution and easier algorithm. 

Fioure 14 illustrates an apparatus based on the Sagnac interferometer including a scann.ng mirror 152 ! and a 
beamsplitter 1 53 In this case, however, the scanning mirror 1 52, rather than the beamsplitter 1 53. is rotated about an 
axis perpendicular to the optical path (and perpendicular to the plane of the paper) to effect the scanning. As shown m 
Fiqure 14. mirror 152 is a double-Sided mirror. 

Figure 15 illustrates a similar arrangement as in Figure 14. except that the scann.ng mirror 162 ,s s.ngle-s.ded and 
is rotated as mirror 152 in Figure 14. to effect the OPD scanning via the beamsplitter 1 63. 

In Figures 14 and. 15 the interferometer is not moving, but it works as in the moving case because of the particular 

arrangement of the scanning mirrors. 

Further variations, modifications and applications of the invention will be apparent. 

It will thus be seen that imaging spectrometers constructed in accordance with the present invention do not merely 
measure the intensity of light coming from every pixel in the field of view, but also measure the spectrum of each pixe 
in a predefined wavelength range. They also better utilize all the radiation emitted by each pixel m the f.eld of view at 
any given time and therefore permit a significant decrease in the frame time and/or a significant .ncrease ,n the sensi- 
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tivity of the spectrometer. Such imaging spectrometers may include other types of interferometers and optical collection 
and focusing systems, and may be used in a wide variety of applications, including medical diagnostic applications, 
remote sensing for geological and agricultural investigations, and the like. 

While the invention has been described with respect to a limited number of embodiments, it will be appreciated that 
many variations, modifications and other applications of the invention may be made. 

Claims 

1 . A method of analyzing an optical image of a scene to determine the spectral intensity of each pixel thereof, com- 
prising: 

(a) collecting incident light from the scene; 

(b) passing the light through an interferometer which outputs modulated light corresponding to a predeter- 
mined set of linear combinations of the spectral intensity of the light emitted from each pixel; 

(c) focusing the light outputted from said interferometer on a detector array; and 

(d) processing the output of the detector array to determine the spectral intensity of each pixel thereof. 

2 The method as in claim 1 wherein said interferometer is of the moving type in which the optical path difference is 
varied to modulate the light by moving the interferometer or an element of the interferometer, such that at each 
instant each detector sees a different point of the scene and its signal is a linear combination of the spectral content 
of the light coming from each pixel, and that when the scanner completes one interferometer scan, the scene will 
have been scanned at all relevant linear combinations of the spectral content and wherein said detector array is 
two-dimensional. 

3 The method as in claim 1 wherein said interferometer is of the moving type in which the optical path difference is 
varied to modulate the light by moving the interferometer or an element of the interferometer, such that at each 
instant each detector sees a different point of the scene and its signal is linear combination of the spectral content 
of the light emitted from each pixel, and that when the scanner completes one interferometer scan and one scene 
scan the scene will have been scanned at all relevant linear combinations of the spectral content, wherein said 
detector array is one-dimensional and further comprising scanning said light in one dimension prior to passing said 
light through said interferometer. 

4 The method as in claim 1 wherein said interferometer is of the non-moving type in which the optical path difference 
varies with the angle of incidence of the incoming light to modulate the light, such that at each instant each detector 
sees a different point of the scene and its signal is a different linear combination of the spectral content of the dif- 
ferent points of the scene, and that when the scanner completes one scan, the scene will have been scanned at all 
relevant linear combinations of the spectral content of the scene, wherein said detector array is two-dimensional 
and further comprising scanning said light in one dimension prior to passing said light through said interferometer. 

5 The method as in claim 1 wherein said interferometer is of the non-moving type in which the optical path difference 
varies with the angle of incidence of the incoming light to modulate the light, such that at each instant each detector 
sees a different point of the scene and a different linear combination of the spectral content of the different points 
of the scene and that when the scanner completes scanning one frame, the complete frame w.ll have been 
scanned at all relevant linear combinations of the spectral content of the scene, wherein said detector array is one- 
dimensional and further comprising scanning said light in two dimensions prior to pass.ng said light through said 
interferometer. 

6. Apparatus for analyzing an optical image of a scene to determine the spectral intensity of each pixel thereof, com- 
prising; 

(a) means for collecting incident light from the scene; 

(b) an interferometer; 

(c) means for passing the light through said interferometer such that it outputs modulated light corresponding 
to a predetermined set of linear combinations of the spectral intensity of the light emitted from each pixel; 

(d) a detector array; 

(e) means for focusing the light outputted from said interferometer on said detector array: and 

(f) means for processing the output of the detector array to determine the spectral intensity of each pixel 
thereof. 
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ThP aonaratus as in claim 6 wherein said interferometer is of the moving type in which the optical path difference 

instant each detector sees a different point of the scene and .ts s.gnal ,s a hnear ~rnb.nat.on of 

of the light coming from each pixel, and that when the scanner comp ietes one ^rferomet ersc aa h ^ c ^ ne w » 

have been scanned at all relevant linear combinations of the spectral content and where.n sa,d detector array 

two-dimensional. 

The apparatus as in claim 6 wherein said interte.oma.er is ot the moving type in «hich *a optical <»* 

^ra^^^ 

light through said interferometer. 
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into said interferometer. 
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(54) Method and apparatus for spectral analysis 

(57) A method of analyzing an optical image of a 
scene to determine the spectral intensity of each pixel of 
the scene, which includes collecting incident light from 
the scene: (b) passing the light through an interferom- 
eter which outputs modulated light corresponding to a 
predetermined set of linear combinations of the spectral 
intensity of the light emitted from each pixel: focusing 
the light outputted from the interferometer on a detector 
array; and processing the output of the detector array to 



determine the spectral intensity of each pixel thereof. If 
the interferometer is of the moving type scanning in one 
dimension is required where the detector array is one 
dimensional, and no scanning when the detector array 
is two-dimensional. If the interferometer is of the non- 
moving type scanning is required in one dimension 
when the detector array is two-dimensional, and in two 
dimensions when the detector array is one-dimensional. 
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